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We present low-temperature anelastic and dielectric spectroscopy measurements on the perovskite 
ionic conductor BaCei_i,Yi,03_a;/2 in the protonated, deuterated and outgassed states. Three 
main relaxation processes are ascribed to proton migration, reorientation about an Y dopant and 
tunneling around a same O atom. An additional relaxation maximum appears only in the dielectric 
spectrum around 60 K, and does not involve H motion, but may be of electronic origin, e.g. small 
polaron hopping. The peak at the lowest temperature, assigned to H tunneling, has been fitted 
with a relaxation rate presenting crossovers from one-phonon transitions, nearly independent of 
temperature, to two-phonon processes, varying as T^, to Arrhenius-like. Substituting H with D 
lowers the overall rate by 8 times. The corresponding peak in the dielectric loss has an intensity 
nearly 40 times smaller than expected from the classical reorientation of the electric dipole associated 
with the OH complex. This fact is discussed in terms of coherent tunneling states of H in a cubic 
and orthorhombically distorted lattice, possibly indicating that only H in the symmetric regions of 
twin boundaries exhibit tunneling, and in terms of reduction of the effective dipole due to lattice 
polarization. 



I. INTRODUCTION 



Perovskite cerates and zirconates are a class of ma- 
terials that, with appropriate doping, exhibit ionic con- 
ductivity both of O vacancies and protons, and therefore 
are suitable as solid electrolytes for fuel cells, gas sen- 
sors and other electrochemical devices. Dissolution of 
H is achieved in two stepsii the material is doped with 
lower valence cations, e.g. BaCei_a;Ya;03_3;/2 (BCY) 
where the partial substitution of Ce'*"'" with intro- 
duces charge compensating O vacancies; the material is 
then exposed to a humid atmosphere at high temper- 
ature, so that the H2O molecules may dissociate, each 
of them filling an O vacancy and contributing with two 
H atoms. Infrared spectroscopy^ and diffraction"^ exper- 
iments indicate that H is bound to an O, forming an 
OH^ ion, but also makes some hydrogen-bonding with 
a next-nearest O atom. The proton diffusion is believed 
to consist of a rapid rotation about the O atom to which 
is associated, and slower jumps to one of the eight next- 
nearest neighbor O atoms with which an instantaneous 
hydrogen bond is established (see Fig. [7^ below). It 
has also been proposed, however, that in some distorted 
zirconates and titanates the rotational barrier may be 
higher than the transfer barrieri^ In various cubic per- 
ovskites ABO3 there are four equilibrium orientations of 
the OH" complex, with an OH separation of 0.9 — 1.0 A 
and H pointing in the (100) directions perpendicular to 
the B-O-B bond'^i^ in the presence of dopants or non 
cubic distortions such positions would be shifted."^ It is 
very likely that the fast local motion of H about the same 
O is dominated by tunneling, but so far no quantitative 
measurements of the associated correlation times have 



appeared, except for quasi-elastic neutron scattering ex- 
periments on hydrated Ba(Cao.39Nbo. 51)02.91, where a 
fast local motion has been detected above room temper- 
ature with an apparent activation energy of ~ 0.1 eV^^ 
Also in SrCeo.95Ybo.05O2.97 a broad quasielastic compo- 
nent has been attributed to fast proton rotation,^ but no 
reliable measurement of the associated rate was possible. 

There is some controversy on the effect of dopants on 
the proton mobility, as reviewed in Ref. [^: on the one 
hand there are several indications for trapping 
with formation of stable dopant-H complexes, but it has 
also been proposed that the excess doped charge dis- 
tributes over all O sites, causing an increase of the hop- 
ping barrier for the proton over the whole latticeJ^ 

Anelastic and dielectric spectroscopies may contribute 
to answer to such issues, since both an electric and elastic 
dipole are associated with a 0H~ ion or Y-OH~ complex, 
and each type of jump or reorientation with characteristic 
time T (T) causes a maximum of the losses at the tem- 
perature T and frequency w/27r such that ujt {T) ~ 1. 
Anelastic relaxation is particularly useful, since it is al- 
most insensitive to electronic conduction and is not af- 
fected by charged interface layers. We present anelastic 
and dielectric spectroscopy measurements on protonated, 
deuterated and outgassed BaCei_2:Y^03_3;/2, where dif- 
ferent relaxation processes are ascribed to H migration, 
reorientation about an Y dopant and tunneling about a 
same O atom; the focus will be on the last type of motion. 

II. EXPERIMENTAL 

The starting powders of BCY were prepared by auto- 
combustion synthesis, which is an easy and convenient 
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solution-based method for the preparation of nanometric 
mixed-oxide powders. The method^** is an improvement 
of that for metal citrates described in the literature.— 
Stoichiometric amounts of highly purified metal nitrates 
were dissolved in distilled water and mixed with citric 
acid which acted both as metal ions complexant and as 
fuel. The citric acid to metal nitrates ratio was main- 
tained to 2, ammonium nitrate was added to regulate 
the fuel to oxidant ratio (citric acid/total nitrate ions ra- 
tio) to 0.4 and ammonia solution (30%wt) was added to 
regulate the pH value at 6. The water solution was left 
to evaporate at 80 °C under constant stirring in a beaker 
immersed in a heated oil bath, until a whitish and sticky 
gel was obtained. The temperature was then raised to 
200 °C until the gel became completely black and dry. 
The beaker was then put directly on the hot-plate at 
250 — 300 °C until the auto-combustion reaction occurred 
leaving the powdered product. Crystallization was com- 
pleted by firing the combusted powders in stagnant air 
at 1000 °C for 5 h. No weight loss occurred during the 
synthesis and no oxide segregation has been detected by 
X-ray diffraction, so that we assume the nominal com- 
positions of BaCei_j;Ya;03_a;/2 with a; = 0.1 (BCYIO) 
and a; = 0.15 (BCY15). The nanopowders were first uni- 
axially pressed at 50 MPa and then isostatically pressed 
at 200 MPa obtaining 60 x 7 x 6 mm bars, which were 
sintered at 1500 °C for 10 h. The bars where cut as thin 
reeds about 4 cm long and 1 mm thick, whose major faces 
where covered with Ag paint. 

The maximum molar concentration Ch/d of H f-nd D 
was measured from the change of weight when the sample 
state was changed between fully outgassed (up to 730 °C 
in vacuum < 10~^ mbar) and hydrated for 1-2 h at 520 °C 
in a static atmosphere of 50 — 100 mbar H2O or D2O 
followed by slow cooling; it was found ch = 0.14 for a; = 
0.15 and ch ^ 0.086 for x = 0.10, slightly less than the 
theoretical maximum ch — x. 

The elastic compliance s (w, T) = s' —is" was measured 
by electrostatically exciting the flexural modes of the bars 
suspended in vacuum on thin thermocouple wires in cor- 
respondence with the nodal lines; the 1st, 3rd and 5th 
modes could be measured, whose frequencies are in the 
ratios 1 : 5.4 : 13.3, the fundamental frequencies of sam- 
ples being a;/27r ~ 2.8 kHz. The elastic energy loss coeffi- 
cient, or the reciprocal of the mechanical quality factor, 

{uj,T) = s" / s' was measured from the decay of the 
free oscillations or from the width of the resonance peak. 
The elastic compliance s is the mechanical analogue of 
the dielectric susceptibility Xj with corresponding 
to tan (5. 

The dielectric permittivity e — e' — ie" was measured 
with a HP 4284A impedance bridge with a four wire 
probe between 3 and 100 kHz in the same cryostat used 
for the anelastic measurements. After depositing the Ag 
electrodes, an intense dielectric relaxation process, iden- 
tified with the motion of charge carriers within a Schot- 
tky barrier at the electrode interface, completely masks 
the true bulk relaxation. Such an effect was suppressed 



by applying 40 V up to ~ 500 K, switching the direction 
of the dc current in order to avoid electromigration of O 
vacancies or protons»ii 

A defect hopping or reorienting with characteristic 
time r and causing a change AA of its elastic quadrupole 
and Ap of its electric dipole contributes with a Debye 
peak to the imaginary parts of the elastic compliance s 
and dielectric permittivity e asi^ 



c (AXf 



LOT 



(UJT) 



and 



c {ApY 



UJT 



3eoWoA:BTcosh^ {E/2T) 1 + {ujTy 



(1) 



(2) 



where c is the molar concentration of relaxing defects, 
Vq the molecular volume, and E is the energy difference 
between the states participating to relaxation. When 
E >T, the higher energy state becomes less populated 
and therefore there is a reduced change of the defect pop- 
ulations under application of the probe field with respect 
to the case E = 0; the consequent reduction of the relax- 
ation strength is described by the factor sech^ {E/2T)J^ 



III. RESULTS 
A. Anelastic spectra 

Figure [T] presents the anelastic spectra of 
BaCeo.gYo.iOs-i measured at the fundamental fre- 
quency of 2.8 kHz: 1) in the as prepared state, therefore 
nearly saturated with H2O; 2) after outgassing H2O at 
900 °C in a flux of pure O2 for 2.5 h; 3) after keeping in a 
static atmosphere of ~ 70 mbar D2O at 520 °C followed 
by cooling at 1 °C/min; 4) after having measured in 
vacuum up to 500 °C with partial loss of D2O. 

There are at least three relaxation processes clearly 
due to the presence of H or D: the peak labeled Py 
at ~ 300 K, Ph at ~ 200 K and Pt at - 30 K; above 
400 K start the contributions due to O vacancies. All 
these peaks shift to higher temperature when measured 
at higher frequency, and a preliminary analysis yields ac- 
tivation energies of 0.58 eV with tq ~ 3 x 10"^'' s for 
Py and ~ 0.4 eV for Ph; the latter is likely composed of 
different processes. These peaks are quite broader than 
pure Debye relaxations and the possible difference in re- 
laxation parameters between H and D is not apparent 
without an accurate analysis; they must be due to hop- 
ping of H between different O atoms. Peak Pt is the focus 
of the present work and shifts from 29 to 38 K when H is 
replaced with D, clearly indicating that it is due to the 
fast motion of H about a same O atom with a dynamics 
dominated by tunneling. An additional peak at ~ 100 K 
might also be due to H and possibly shifts to higher tem- 
perature after substitution of H with D, but its nature 
is not as clear as for the other processes and we will not 
discuss it further. 
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FIG. 1: Anelastic spectra of BCYIO measured in different 
conditions of hydration with H2O and D2O. 




FIG. 2: Dielectric (left-hand scale) and anelastic (right-hand 
scale) spectra of BCY: e" of hydrated BCY15 at 3, 10 and 
30 kHz (continuous lines) and e" of hydrated BCYIO at 3 kHz 
(curve 1). The dotted lines are Q^^ of hydrated (2) and 
outgassed (3) BCYIO measured at 2.8 kHz. 



B. Comparison between anelastic and dielectric 
spectra 

All the peaks appearing in the anelastic spectrum are 
present also in the dielectric one, which also displays an 
additional maximum Pd at ~ 60 K. Figure [5] shows e" 
of hydrated BCY15 and BCYIO (ch = 0.078, curve 1); 
the latter is compared with the elastic counterpart 
in the hydrated (curve 2) and outgassed (curve 3) states, 
measured at the same frequency of 3 kHz. It appears 
that the intensities of peaks Pt, Ph and Py span more 




FIG. 3: Low-temperature region of the dielectric and anelastic 
spectra of BCY of Fig. [1 The e" (left) and (right) 
scales are chosen in order to emphasize the near coincidence 
of dielectric and anelastic Pt (curves 1 and 2). 



than four orders of magnitude in e" but only two in ^ 
meaning that the fast tunneling motion of H produces a 
much smaller change in the electric dipole than in the 
elastic dipole, compared to the hopping motion. On the 
other hand, the peak temperatures and activation en- 
ergies appear to be practically the same in the dielec- 
tric and anelastic losses. These dielectric measurements 
on BCY15 have been made after full elimination of the 
charge relaxation at the electrodes, so that the intensi- 
ties of Py and Ph are reliable. Figure [3] shows e" of 
both hydrated (thick lines) and deuterated (thin lines) 
BCY15. While peak Pt shifts to higher temperature after 
the isotope substitution of H with D as for the anelastic 
case, peak Pd remains completely unaffected; this fact, 
together with its absence in the anelastic losses, indicate 
that Pd is not connected with the H motion, but is likely 
of electronic origin. Curves 1—3 in Fig. [3]are the same di- 
electric and anelastic curves of BCYIO appearing in Fig. 
[21 but this time the two scales are chosen in order to yield 
the same intensity of Pt. Unfortunately, experimental 
difficulties connected with the elimination of the metal- 
semiconductor barrier at the electrodes did not allow us 
to obtain significant dielectric spectra in the outgassed 
state, so that it is possible to determine that the inten- 
sity of Pd almost triplicates in the hydrated state passing 
from 10% to 15%Y doping, but nothing can be said on 
its dependence on the content of H or O vacancies. 
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IV. DISCUSSION 
A. Proton hopping and trapping 

In what follows we will identify peak Py with hopping 
of H among the O atoms of YOg octahedra, namely with 
the reorientation of the Y-H complex, and Ph with hop- 
ping over CeOg octahedra, with the possible contribution 
of the formation/dissociation of Y-H complexes. These 
assignments are suggested by the fact that it is natu- 
ral to assume that (0H)+ ions, having a formal charge 
+1 with respect to 0^~, may form relatively stable com- 
plexes with trivalent dopants Y"^"*", having formal charge 
-1 with respect to Ce*+. Then, Py with higher intensity 
and activation energy should be due to the more numer- 
ous H atoms associated with Y, while Ph should be due 
to the faster jumps of the less numerous H atoms not as- 
sociated with Y. The assignment of anelastic relaxation 
peaks to dopant-H complexes has already been proposed 
for BaCei_^Nd:,03-i»^ and (Ba,Sr)Cei_a;Yb:,03, - al- 
though the argument that the (0H)~ ion has the same 
symmetry as the crystal and therefore cannot produce 
anelastic relaxation unless forming defect complexes^S is 
incorrect. Additional experimental indications that H is 
trapped by trivalent dopants in BaCei-ajY^-Os are the 
analysis in terms of two components of the quasielastic 
neutron scattering peak,^*' and the observation with EX- 
AFS of an enhancement of the disorder in the environ- 
ment of Y after hydration.— Also first-principle calcula- 
tions and Monte-Carlo simulations of the proton diffusion 
indicate that dopants act as trapsi^i^ 

On the other hand, there are also experiments sug- 
gesting the absence of significant trapping, like an NMR 
investigation— on BaCei_a;Ya;03 with x = 0.01 and 0.1, 
where the correlation time deduced from the "'^H NMR re- 
laxation is the same at both doping levels and reproduces 
the conductivity with a simple hopping model without 
trapping. At this stage, it cannot be completely excluded 
from our data that there is indeed very little trapping 
effect from Y dopants, and the two main anelastic re- 
laxation processes Py and Ph are associated with H hop- 
ping among 01 and 02 atoms of the orthorhombic struc- 
ture, having different symmetries. It is possible that H 
binds to 01 and to 02 with different probabilities, as 
neutron diffraction experiments and molecular dynamics 
simulations suggest)^i23, and jumps within the respective 
sublattices with different rates, so giving rise to peaks 
with distinct intensities and temperatures. In the dis- 
cussion we will also mention the possibility that in the 
distorted orthorhombic structure the reorientation of the 
OH ions is much slower than found in the higher temper- 
ature phases, and Ph is due to such slow reorientation. 



B. Fit of the anelastic Pt 

The relaxation modes of the fast reorientation of OH 
among four positions can be found by solving the rate 



equations for classical hopping between the four H sites 
or the quantum mechanical problem with tunneling be- 
tween nearest neighbor orientations. In the classical sym- 
metric case there are two modes: one active in the dielec- 
tric relaxation and another in the anelastic relaxation 
with a rate twice larger; in the case that the site ener- 
gies are different, however, the modes become three, all 
contributing to both anelastic and dielectric relaxation, 
and with possibly widely changing rates, depending on 
the type of asymmetry and the degree of coherence of 
the eigenstates of H plus polaron-like distortion of the 
surrounding atoms. At variance with the geometrically 
analogous case of the reorientation of Zr-H complexes in 
Nb^i^i2£i22 in highly doped and orthorhombic BCY it 
is not possible to distinguish different peaks arising from 
these modes, and therefore we will limit ourselves to fit 
Pt with a single relaxation time r plus broadening, as 
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where the parameters a, /3 < 1 produce broadening of the 
low- and high-temperature sides of the peak, respectively; 
when a = (5 = 1 the above expression reduces to a Debye 
peak. The parameter E is the energy difference between 
the configurations involved in relaxation, and must be 
introduced in order to reproduce the enhancement of the 
peak height at higher frequency. The relaxation rate 
was modeled as 

= T-pi^coth(£;/2fcBT) + (T/r2ph)"+ (4) 
+ro-i exp yW/T) . (5) 

Such an expression is similar to that used by Kuskowsky, 
Lim and Nowick)^ for the low temperature dielectric 
relaxation in Bai_a;Nd2;Ce03_a;/2; it does not rely on 
a model of the interaction between a precise defect ge- 
ometry and the actual phonon bath, but is able to de- 
scribe the main hopping regimes, including tunneling 
in an insulating crystal. One expects that, starting 
from low temperature, the transitions between the de- 
fect eigenstates occur through processes involving one- 
phonon, then two-phonons and finally several phonons 
or semiclassical hoppingi^ In the one-phonon regime 
the transition rate is coth {E/2kBT), which becomes 
temperature independent when fceT" is smaller than the 
separation E between the eigenstates; the latter has a 
form of the type E ~ + where the tunneling ma- 
trix element t is expected to be smaller than the typi- 
cal energy asymmetry a between the site energies due to 
the orthorhombic distortion and disordered nature of the 
BCY solid solution; in the fit i? = a. The two-phonon 
relaxation rate approximately depends on temperature 



through a power law 



29.30.31 



with 5 < n < 9 depend- 



ing on the type of interaction with acoustic phonons and 
on the energy difference between H sites. For interac- 
tion with optical phonons of frequency wq an Arrhenius- 
like exp {—huo/kBT) dependence is predicted;^ At suf- 
ficiently high temperatures an Arrhenius-like temper- 
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FIG. 4: Fit of the anelastic Pt of hydrated and deuterated 
BCYIO, measured at 2.9 and 15 kHz. 
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FIG. 5: Relaxation rates used in the fits of Fig. U The 
rate for deuterium is decomposed into one- , two-phonon and 
Arrhenius-like contributions. 



ature dependence is always found, although the pre- 
exponential term tq is not necessarily related to the fre- 
quency of a local H vibration mode, but depends also on 
the overlap of the wave functions of H in adjacent sites, 
and the barrier W may contain significant corrections 
due to plionon fluctuations with respect to the static po- 
tential. 

We found that all the three contributions to are 
necessary to obtain a good fit: the one-phonon term re- 
produces the low temperature broadening of Pt without 
resorting to extremely small values of a. The introduc- 
tion of the power low significantly improves the fit quality 



around the maximum. The Arrhenius contribution is not 
essential to obtain good fits, since its suppression can be 
partially compensated by an increase of /3, namely by di- 
minishing the broadening at high temperature; however, 
if one wants to keep a ~ /3 the Arrhenius contribution 
must be included. The condition a ~ /3 appears desir- 
able, if such parameters should describe the broadening 
of a Debye peak due to lattice disorder, and not to other 
types of interactions, like collective interactions among H 
atoms. 

The continuous lines in Fig. 3] are fits of anelastic Pt 
of BCYIO both hydrated and deuterated with the above 
expressions plus a linear background. In view of the par- 
tial duplication of the effects of some parameters, like 
the high temperature broadening /3 and the Arrhenius 
To and W, these fits are not unique, and we tried to 
obtain a physically sound combination of the parame- 
ters. The mean asymmetry energy is determined quite 
precisely as E'/Zcb = 64 K from the temperature depen- 
dence of the peak intensities. Regarding broadening, it 
is possible to obtain good fits with a = f3 = 0.5 for both 
H and D, although (3 = 0.38 for H gives a slightly bet- 
ter interpolation, as in Fig. [U these values of a and 
f3 imply a broadening that is conspicuous but expected, 
in view of the high lattice disorder. The other parame- 
ters are: To = 7.7 x IQ-i" s (1.1 x IQ-^^)^ W = 7U K 
(820), T2ph = 6.7 K (8.4), n ^ 7, t'J^ = 7300 s^^ 
(990) for H (D). The resuhing t^\^ (T) are plotted in 
Fig. [5] and it turns out that Tjj^/t^^ ~ 8 at all tem- 
peratures; such a ratio certainly indicates a non-classical 
effect of the isotope mass on the H dynamics. The bar- 
rier W ~ 800 K of the Arrhenius-like contribution is close 
to the activation energy for rotational diffusion obtained 
from quantum molecular dynamics simulations^i^^ and 
to that extracted from quasi-elastic neutron scattering 
in Ba(Cao.39Nbo.6i)02.9i'^ It cannot be excluded, how- 
ever, that it rather originates from two-phonon interac- 
tion with optical phononSf^ considering that the infrared 
absorption bands in various cerates range from hw /k-Q = 
600 K to 1000 K (400 to 700 cm-^)^ The power law 
with n = 7 instead is typical of two-phonon transitions 
of asymmetric states interacting with acoustic phonons; 
lower values of n yield definitely worse fits, while n ~ 7.8 
is found if the Arrhenius contribution is omitted. 



C. Fit of dielectric Pt and Pd 

The most reliable fit is on BCYIO, where peak Pd has 
a reduced intensity. Below 100 K, e' {us, T) has already 
approached the limiting high frequency value ^oo = 20.1, 
so that the e" (T) = x tan 5 (T) and in Fig. [HI we 
show tan (5 at 10, 30 and 100 kHz; at lower frequency 
the noise was too large to add significant information. 
Note that, due the imperfect compensation of the cables 
impedance at low temperature, and the very small values 
of the dielectric losses, the tan 5 curves can be arbitrarily 
shifted in the ordinate scale. In fact, it was verified by 
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and Td = Tdo exp (Wd/T). The thick curves in Fig. [6] are 
the resulting fit with hnear backgrounds, and also Pt and 
Pd are shown. The parameters of Pd are: ad = 150 K 
(150), Tdo = 1 X 10-15 s (6 X 10-15), Wd = 1270 K (1270) 
and Qfd — 0.26 (0.38), where the values in parenthesis 
are obtained from the measurements on BCY15 (data of 
Fig. The small values of Tdo are in agreement with 
the electronic origin of the relaxation, but, considering 
the extreme peak broadening (small ad), it cannot be 
excluded that a correlated dynamics is present, which 
may be better described by a Vogel-Fulcher type t (T) 
rather than Arrhcnius. 



D. The intensity of Pt 

1. The electric dipole of the OH group 



FIG. 6: Fit of the dielectric spectrum of BCYIO. 



switching on and off the cable compensation that below 
150 K such a correction became independent of temper- 
ature and therefore introduced only a shift of the tan 6 
curves. 

The shapes and temperatures of dielectric and anelas- 
tic Pt need not to be identical, since the various re- 
laxation modes contributing to them may have different 
strengths for anelastic and dielectric relaxation. Yet, we 
chose not to let all the parameters of Pt to vary freely, due 
to the still overwhelming presence of Pd at higher temper- 
ature and possibly to the presence of another relaxation 
at lower temperature. Therefore, all the parameters of Pt 
where set to the values of the anelastic fit of Fig. 01 with 
a factor multiplying the anelastic r^^ as the only degree 
of freedom: t^^j = r r^^; it turns out that r — 10.6 
yields a good fit. Peak Pd, appearing only in the di- 
electric spectrum, has an intensity Aed that clearly rises 
with temperature (see Fig. [3]) and there are two main 
possible causes for such a behavior: i) Pd is due to the 
relaxation of a thermally excited state with energy E over 
the ground state, e.g. charges from a ionized defect, and 
therefore the relaxation strength contains the population 
of the ionized state as a factor, Aed oc 1/ cosh(i?/2T); ii) 
relaxation occurs between two states differing in energy 
by a, so that Asd oc nin2 = 1/ cosh^ (a/2T). It is also 
possible that both mechanisms are present, but the data 
do not allow to distinguish between these possibilities, 
and we will consider the case a 0, E = 0; note that the 
two mechanisms give similar temperature dependence of 
Asd in a broad temperature range if E ^ 1.5a. The 
shape of Pd could be better reproduced with the Cole- 
Cole expression for broadening, so that the expression for 
fitting Pd was 

„ _ Ad sin (7rQ;d/2) 

'^^ T sech^ (ad/2T) cosh [ad In (wTd)] + cos (7rad/2) 



The dependence of the intensity of the anelastic Pt on 
the H content and the marked shift to higher tempera- 
ture with the heavier D isotope mass leave little doubt 
that anelastic Pt is due to the fast motion of H around a 
same O with a dynamics dominated by tunneling. It is 
also clear that Pt has its dielectric counterpart (Fig. [31 
curves 1 and 2). A puzzling feature of the dielectric Pt 
is its very small intensity, and therefore we discuss now 
about the estimated strengths of dielectric and anelastic 
relaxation. The dielectric relaxation strength associated 
with various types of H jumps should be relatively easy to 
estimate in a medium which is not particularly highly po- 
larizable like BCY. Let us first consider Py, interpreted 
as the reorientation of an effective dipole Y'— OH* where 
the Kroger- Vink notation expresses the fact that Y'^+ 
has an excess — e charge with respect to Ce^"*" of the per- 
fect lattice and OH- has a +e charge with respect to 
O^^. On the time scale of the H reorientation among 
different faces of the cube containing Y, the fast motion 
around each O is completely averaged out with barycen- 
ter near O, so giving rise to an effective dipole — ea/2, 
where a is the lattice constant of the cubic perovskite 
(Fig. [7^). A jump to a different cube face will cause a 
change of the dielectric dipole by Ap = eV^a; setting^ 
a = 4.39 A, Eq. ((2]) gives a dielectric relaxation strength 
Ae^ — cy-h X 1060 at 300 K; therefore the intensity 
Ae ~ 30 of Py at X = 0.15 (Fig. HD is obtained setting 
cy-h — 0.03, namely that at this temperature ~ 20% of 
H is trapped by Y. It is likely, however, that cy-h is 
closer to ch but the effective dipole is reduced by the 
surrounding polarization and distortion. The dipole in- 
volved in Pt instead, is associated with the OH ion in its 
four equilibrium positions (Fig. Uh)- 

P^"" = -P?" = 9rfx, (6) 

p?" = -Pr = -<idy, 

where c? ~ 1 A is the 0-H separation. The actual charge 
on O and H is likely reduced by polarization effects, 
and it has been calculated by atomistic simulations^ 
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FIG. 7: (a) Network of the H sites in the cubic structure 
of BCY; Ba — white, O = gray, H sites = black, Y/Ce at 
the centers of the octahedra. (b) Electric dipoles pp^ and 
elastic quadrupoles Ap^ associated with the four orientations 
of the 0H~ ion . (c) Distorted environment of a set of four 
sites around a 02 atom in the orthorhombic structure (atomic 
positions except H from Ref.— in Pnma setting). 



as — 1.426e and 0.426e, so giving q = 0.426e. A reori- 
entation of the OH ion by 90° causes Ap'^^ — q^/2d 
and, setting Ea — 64: K from the fit, the resulting re- 
laxation strength of peak Pt at T ~ 35 K should be 
Ae^ = cy-h X 41. It results that Pt at 35 K should 
be ~ 26 times smaller than Py at 300 K, but instead it 
is nearly 1000 times smaller, therefore nearly 40 times 
smaller than expected. 



2. The elastic quadrupole of the OH group 



of the strain associated with the OH ion and cannot be 
estimated in an obvious manner as the electric dipole. It 
can only be said that the OH ion in the cubic configura- 
tion of Fig. [7] has orthorhombic symmetry and therefore 
in a xy cube face it is 



\OH 



vOH 















\OH 
^^3 



(7) 



and A^J have the x and y components exchanged. The 
reorientation of OH by 90° causes a change AA*-*^ = 
||Af^ — A^^||. Onthe longer time scale of relaxation Py , 
the four positions of the face _L z produce a tetragonal 
quadrupole 



i (A?H + A?H) 

A^ = I i (A?H + AOH) 

AOH 



(8) 



plus, in case of Y-OH complex, a possible additional dis- 
tortion 




(9) 



Then, the anelastic peak Py associated with the reori- 
entation of the tetragonal defect complex Y-OH, has 
AAY = (AY - Xj) + A™ _ 1 (aOh + aOH) . Lacking any 
method for estimating these components of the elastic 
quadrupoles associated with OH and Y-OH complexes, 
we cannot exclude that AA*-*^ is ~ 10 times smaller than 
AA^, so determining a reduction of the intensity of Pt by 
two orders of magnitude with respect to Py. Then, the 
small intensity of the anelastic Pt does not necessarily 
constitute a problem as the dielectric intensity does. 



Anelastic and dielectric relaxation strengths of tunneling 
states in a cubic environment 



The explanation why the intensity of the dielectric Pt 
is so much smaller than expected from a simple estimate 
of the magnitude of dipole change would be easy if BCY 
was cubic or very close to cubic. In fact, in case of coher- 
ent tunneling among four nearly equivalent positions, the 
effective dipole strength for transitions between H eigen- 
states would be much smaller than for hopping among 
the same positions, at variance with the anelastic case. 
The Hamiltonian of the symmetric four-level tunnel sys- 
tem (FLS) of H may be written in the basis \i) where H 
is localized in each site i = a — d as: 



H, 



cub 



t t 

t t 

t t 

t t 



(10) 



The anelastic case is different, since the elastic 
quadrupole A corresponds to the long range component 



where t is the effective tunneling matrix element between 
adjacent sites dressed with the interaction with phonons 



8 



and the site energies are all set to 0. The eigenstates of 
i?cub are 



|1,4) 



" ±1 ■ 




" ■ 




" -1 ■ 


1 

±1 


, |2)(x 


-1 



, |3)« 



1 


1 




1 








(11) 



two states are delocalized over all sites, with energies 
i?i,4 = =Fi/2j and two are delocalized over either pair of 
opposite sites, with energies i?2,3 — 0. The associated 
electric dipoles p^'"^ and elastic quadrupoles A^^^ may 
be obtained taking the matrix elements Xf^^ = {i\ A*^^ \i) 
and pf^^ = {i\p'^^ with A*^^ and p*-*^ given above. 
It is evident that p'-'^ averaged over any of the above 
states is null, since opposite sites are equally occupied 
and cancel out the respective dipoles. Then, no dielec- 
tric relaxation is expected from H tunneling in a cubic 
environment, or also next to a Y dopant, which leaves the 
fourfold symmetry of the FLS. The elastic quadrupoles, 
instead, being centrosymmetric as in Eq. ([7]), are equal 
within pairs of opposite sites, and transitions between 
the two intermediate eigenstates cause a change of elas- 
tic quadrupoleby AA°^ = ||A°^ — Aj*^!]; transitions be- 
tween states 2,3 and 1,4 cause a change by AA*-'^/2 and 
those between 1 and 4 no anelastic relaxation. There- 
fore, the formation of FLS with weak deviations from 
fourfold symmetry would explain why the dielectric re- 
laxation from H tunneling is much more effectively sup- 
pressed than the anelastic one, with respect to classical 
reorientation. 

A picture like this, with H and D performing coherent 
tunneling within nearly symmetric FLS, has been thor- 
oughly studied by anelastic relaxation in Nb with sub- 
stitutional trapsi2^i2^i2^i21 Also in that case H tunnels 
within rings of four equivalent tetrahedral sites on the 
faces of the bcc cells with a dopant in the center, and per- 
forms overbarrier jumps to the neighboring rings, which 
form a network exactly as in the cubic perovskite (see Fig. 
[7^). In Nbi„2:Zr2; single crystals with x ~ 0.0013 it was 
also possible to distinguish well separated relaxation pro- 
cesses arising from transitions involving eigenstates with 
different symmetries, only moderately perturbed by in- 
teractions among dopants?^ in addition, the dependence 
of the anelastic relaxation due to the slower reorienta- 
tion among different cube faces on the symmetry of the 
excitation stress^ provides evidence that the symmetry 
of the FLS persists at least up to 150 K, implying that 
the H eigenstates maintain coherence up to that tem- 
perature. It should be noted that the measured effective 
tunneling matrix element between tetrahedral sites in Nb 
is t 0.2 meV for H and 0.02 meV for D, corresponding 
to t/kB = 2 and 0.2 K or t//i = 4 x 10^° and 4 x 10^ s^^ 
respectively; it appears therefore that coherence is main- 
tained both at temperatures orders of magnitude larger 
than i/fce and with tunneling frequencies of the order of 
phonon frequencies. 

A theoretical basis for the persistence of quantum 
coherence at such high temperatures comes from ana- 



lytical and numerical analysis of the centrosymmetric 
FLS, where coherent oscillations of the H populations 
are found even for strong interaction with the thermal 
bath,"^^ and from the analysis of the dynamics of polarons 
using the dynamical mean-field approximation,"^^ which 
is a non-perturbative approach. It is shown that the co- 
herence of the state including tunneling particle and sur- 
rounding polaronic distortion is maintained to temper- 
atures up to a substantial fraction of the energy of the 
phonon coupled with the particle, k^T ~ 0.2fia;o, also 
in the case of strong coupling, where the energy for the 
polaron formation is Ep > HajQ. We assume i?p ~ 1 eV 
from the theoretical estimate^ of the self-trapping en- 
ergy of H in BaZrOs, and that the strongest coupling 
is with the 0-Ce-O bending mode (in BaCeOa hujQ = 
41 meV"^^) modulating the distance with the neighboring 
O atoms and therefore the hydrogen bonds with them. 
Then we are in the limit of strong coupling, and the dy- 
namical mean-field analysis"^^ ensures us that coherent 
states may be maintained up to 0.2fiUJo/kB — 100 K; this 
holds for both H and D, since the smallness of the dressed 
tunneling matrix element should not be a problem, un- 
til the tunneling frequency is larger than the measuring 
frequency. 



4- Tunneling states in the orthorhombic lattice 

The main problem with the above explanation of the 
smallness of the dielectric relaxation strength of Pt is that 
BCY at low temperature is not cubic but orthorhom- 
bic, and the FLS should be far from symmetric. Figure 
[Tj; shows an undistorted ring, as in the cubic case, be- 
tween two octahedra tilted as in the orthorhombic Pnma 
low temperature structure of BaCeo.gYo.iOa-a."^ The O 
atoms split in two types: 01 at the vertices of the octa- 
hedra along the b axis and 02 near the ac plane; neutron 
diffraction indicates that, at 4.2 K, H in BCY occupies 
a site near the one labeled as d in Fig. [T):, which is also 
the one with the largest distances from the next nearest 
neighbor O atoms. It is therefore reasonable to assume 
that in the distorted orthorhombic structure H occupies 
sites slightly displaced from those in the cubic cell, and 
that the site energy is mainly determined by the mean 
distance I from the two next nearest neighbor O atoms, 
with which some hydrogen bonding can take place. ^ For 
the sites labeled a — d in Fig. [71d it is 7 = 2.66, 2.71, 
2.53, 2.50 respectively, so that the site energies may be 
written as Ea — Eb = a/2, Ec — Ed = —a/2, and the 
Hamiltonian of the tunnel system becomes: 



-f^ortho — 2 



a t 

t a 

t 

t 



t 

t 

—a t 

t —a 



(12) 
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with eigenstates 





a — S 




" a + S ' 


1,2) cx 


Tia-5) 
t 


, |3,4)cx 


T{a + S) 
t 



and energies £'1^2 = {—5 T t) /2, i?3,4 — [5 ^t) /2, where 
5 = \/ a? + i^. This means that, when a 3> t, there 
are two low energy eigenstates with H mainly delocal- 
ized over sites c and d and two higher energy eigenstates 
delocalized over sites a and 6; these states have an electric 
dipole p ~ ed/ oriented roughly midway between the 
two occupied sites. Therefore, while in the limit a <^ t , 
valid for the cubic ideal case, the H atom is delocalized 
over all sites and the electric dipole is averaged out to al- 
most zero, H in BCY at low temperature is expected to 
be in the opposite limit t where the averaging effect 
of the electric dipole occurs only within the pairs ab and 
cd of low- and high-energy states, without a suppression 
of the dipole magnitude. 

The present data do not allow an estimate of t, neither 
do we know of estimates of the energy asymmetry a due 
to the orthorhombic distortion, but based on the com- 
parison with the better known case of Zr-H complexes in 
Nb we should be in the limit a ^ t. In fact, t should 
be smaller than in Nb, since, assuming an 0-H bond 
~ 0.93 A long, the distance between neighboring sites 
is ~ 1.3 A in BCY while it is only 1 A in Nb; in addi- 
tion, the maximum of the anelastic relaxation is shifted 
to higher temperature with respect to Nbi_a;Zr2,Hj,, in- 
dicating slower transition rates; therefore it should be 
t/kB <C 1 K, the value found in Nb. On the other hand, 
we expect a/ fee ^ 10^ K, considering that the random 
strains due to < 1 at% of impurities in Nb cause a/fce 
of tens of kelvin,^^ and the strain associated with the oc- 
tahedral tilts in orthorhombic BCY is certainly larger. 
Another indication that a due to the orthorhombic dis- 
tortion is large comes from the theoretical estimate'^^ 
Ep ~1 eV of the self-trapping energy of H in BaZrOa. 

5. Symmetric tunneling states within twin walls 

After these arguments, it is puzzling that the inten- 
sity of the dielectric Pd in orthorhombic BCY is so 
small. The situation is similar to that found in hydrated 
BaCei_a;Nda;03_5, where a low temperature dielectric 
relaxation exists whose rate has a temperature depen- 
dence indicating tunneling; however, it was concluded 
that if all the protons were responsible for such a re- 
laxation, the intensity should be 50 times larger, and 
therefore only special defect configurations contributed 
to that relaxation.^'* In the present case, unless the OH 
dipole is smaller than estimatedySi less than 3% of the H 
atoms should contribute to Pt and it should be explained 
what kind of particular configuration exhibits tunneling 
and what is the dynamics of the majority H atoms. In 
fact, both neutron spectroscopy at high temperature?'* 



and simulations^ indicate that the reorientation of the 
OH ion is much faster than the hopping between differ- 
ent O atoms in perovskite cerates. A possible scenario 
is that the fast reorientation occurs only in high tem- 
perature phases that are cubic or less distorted than the 
orthorhombic phase, whereas in the latter H is nearly lo- 
calized at lowest energy site, close to site d in Fig. \7i). 
Then the reorientation of the OH ion would be slower 
than that producing Pt and might be identified with Ph 
or some broad peak masked by the tail of Py and by 
Ph, perhaps the one around 100 K in Fig. [TJ In this 
scenario the tunneling motion would appear only in par- 
ticularly symmetric environments, 6.17. at the boundaries 
between different structural domains. We are not aware 
of any study of the density, width and morphology of 
the domain boundaries in BCY, but indirect support to 
this mechanism comes from a simulation on orthorhom- 
bic CaTiOs, where the twin walls are found to be about 
6 pseudocubic cells wide and to trap the O vacancies;'^^ it 
is also proposed that the diffusion of O vacancies should 
be faster within the twin planes, which are more sym- 
metric than the orthorhombic bulk. Also in the case of 
H in BCY, the reorientation rate of H might be faster at 
the twin walls, but the relevance of this effect to the long 
range mobility would be limited, since the rate limiting 
step is not the OH reorientation but the hopping to a 
different O atom. We think however that the effect of 
the greater symmetry at the twin walls should be stud- 
ied also in relation with the hopping mechanism and in 
the high temperature phases of BCY. Also simulations on 
BaZrOs and CaTiOa suggest that the octahedral tilting 
is essential in determining H site energies and diffusion 
paths. ^ 

If H tunneling in the orthorhombic phase indeed occurs 
only in the twin boundaries, than the intensity of peak 
Pt would depend on their density, which in turn may de- 
pend on microstructure and thermal history. This might 
be the reason why the dielectric Pt in BCY15 is nearly 
three times more intense than in BCYIO, a feature oth- 
erwise difficult to explain; additional measurements are 
necessary to clarify this issue. 

Another factor that may contribute to lower the dielec- 
tric relaxation strength of Pt is the lattice polarization 
around the OH^ ion; in fact, the actual dipole is not the 
bare 0H~ dipole with nominal charges ±1, but it is due 
to the 0H~ complex plus the shifted surrounding atoms 
and with charge transfers with respect to the purely ionic 
case. We estimated the intensity of dielectric Pt assum- 
ing the OH dipole calculated in Ref?2i, but the actual 
dipole may be even smaller. 



V. CONCLUSIONS 

Three main relaxation processes have been identified 
both in the anelastic and dielectric spectra as due to hop- 
ping of H around an Y dopant (peak Py), hopping far 
from dopants (peak Ph) and tunneling within the four 
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sites around a same O atom (peak Pt). An additional di- 
electric relaxation maximum around 60 K (Pd) does not 
involve H motion, but rather appears as a relaxation of 
electronic origin like small polaron hopping. Peak Pt can 
be fitted assuming a relaxation rate that is Arrhcnius- 
like above ~ 50 K, possibly due to two-phonon tran- 
sitions with optical phonons rather than to overbarrier 
hopping, below 50 K exhibits a dependence typical 
of two-phonon transitions with acoustic phonons, and fi- 
nally becomes nearly constant below ~ 20 K, as expected 
from one-phonon transitions. The isotopic substitution 
with D decreases the rate by a factor of 8. The dielectric 
spectra are more difficult to analyze, due to the presence 
of peak Pd and possibly other peaks at lower temper- 
ature, but reasonable fits are obtained using the same 
anelastic expression of Pt with a rate increased by ~ 10 
times; such a difference in the rate may be due to the 
fact that, although not clearly distinguishable, there are 
at least three anelastic and dielectric relaxation modes 
contributing to Pt, having different strengths and rates. 

The intensity of the dielectric Pt is nearly 40 times 
smaller than expected from a simple estimate of the re- 
orientation of the electric dipole associated with the OH" 
ion. It is shown that, while this suppression of the di- 



electric relaxation would be easily explained in terms of 
coherent tunneling of H around O in a cubic environ- 
ment, the same seems not to be true in the presence of 
the low-temperature orthorhombic distortion. As alter- 
native or concomitant explanations, it is proposed that 
H tunneling may occur only in the more symmetric cells 
within twin walls, while slower scmiclassical reorientation 
of OH" would occur within the orthorhombic domains; 
in addition, the total dipole of OH" ion and surrounding 
lattice polarization may be smaller than expected. 
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